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Simulating  fire  whirls 

Francine  Battaglia  f,  Kevin  B McGrattan,  Ronald  G Rehm  and  Howard 
R Baum 

National  Institute  of  Standards  and  Technology,  Building  and  Fire  Research 
Laboratory,  100  Bureau  Dr  Stop  8640,  Gaithersburg,  MD  20899-8640 

Abstract.  A numerical  investigation  of  swirling  fire  plumes  is  pursued  to  understand 
how  swirl  alters  the  plume  dynamics  and  combustion.  One  example  is  the  ‘fire  whirl’ 
which  is  known  to  arise  naturally  during  forest  fires.  This  buoyancy-driven  fire  plume 
entrains  ambient  fluid  as  heated  gases  rise.  Vorticity  associated  with  a mechanism 
such  as  wind  shear  can  be  concentrated  by  the  fire,  creating  a vortex  core  along  the 
axis  of  the  plume.  The  result  is  a whirling  fire.  The  current  approach  considers  the 
relationship  between  buoyancy  and  swirl  using  a configuration  based  on  fixing  the  heat 
release  rate  of  the  fire  and  imposing  circulation.  Large  eddy  methodologies  are  used 
in  the  numerical  analyses.  Results  indicate  that  the  structure  of  the  fire  plume  is 
significantly  altered  when  angular  momentum  is  imparted  to  the  ambient  fluid.  The 
vertical  acceleration  induced  by  buoyancy  generates  strain  fields  which  stretch  out  the 
flames  as  they  wrap  around  the  nominal  plume  centerline.  The  whirling  fire  constricts 
radially  and  stretches  the  plume  vertically. 


1.  Introduction 

Fire  whirls  are  a rare  but  potentially  catastrophic  form  of  fire.  These  swirling  buoyant 
fire  plumes  are  known  to  increase  the  danger  of  naturally  occurring  or  post-disaster 
fires.  Researchers  have  spent  the  last  four  decades  investigating  these  devastating  fires 
[1,2, 4-8].  In  order  for  a fire  whirl  to  exist,  there  must  be  an  organized  source  of  angular 
momentum  to  produce  the  large  swirl  velocities  as  air  is  entrained  into  the  fire  plume 
[1,  2],  These  vorticity-driven  fires  occur  over  a large  range  of  length  and  velocity  scales, 
and  significantly  alter  the  entrainment  and  combustion  dynamics. 

One  of  the  earliest  studies  to  mathematically  quantify  a fire  whirl  from  experiments 
was  conducted  by  Emmons  and  Ying  [1].  They  reproduced  a fire  whirl  from  a pool 
of  acetone  centered  within  a rotating  cylindrical  screen.  The  flame  was  observed  to 
lengthen  vertically  and  tighten  radially  with  increasing  screen  rotation.  The  physical 
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explanation  was  that  the  vertical  acceleration  induced  by  buoyancy  generates  strain 
fields  which  stretch  out  the  flames  as  they  wrap  around  the  nominal  plume  centerline. 
A review  by  Morton  [2]  describes  the  underlying  mechanisms  which  produce  and  sustain 
a fire  whirl.  A mathematical  description  of  a swirling  buoyant  plume  arising  from  a point 
source  of  heat  was  approached  by  Thomas  and  Takhar  [3].  Their  analyses  included  cases 
for  which  circulation  decays  to  zero  outside  of  the  plume  or  remains  a finite  constant. 
While  the  analysis  of  these  cases  showed  qualitatively  different  behavoir,  neither  would 
explain  the  observed  plume  dynamics  found  by  Emmons  and  Ying. 

The  extent  of  full-scale  experimental  investigations  of  fires,  in  general,  are  limited 
by  concerns  for  safety,  economy,  control  and  feasibility.  Scale  modeling  is  one  approach 
towards  investigating  fires  under  more  manageable  conditions  by  conducting  small-scale 
experiments  [4,  5,  6].  The  literature  includes  reports  by  Williams  [4]  which  address 
the  importance  of  scaling  mass  fires.  Soma  and  Saito  [5]  have  been  successful  with 
developing  scaling  laws  to  model  prototype  fire  whirls.  An  extensive  review  on  the 
progress  of  experimental  investigations  in  fire  spread  scenarios  has  been  reported  by 
Hirano  and  Saito  [6]. 

Recently,  an  experimental  study  of  swirling  fires  was  pursued  by  Satoh  and  Yang  [7]. 
They  followed  these  experiments  with  a numerical  investigation  as  an  additional  means 
to  substantiate  their  observations  [8].  The  fire  took  place  in  a rectangular  compartment 
with  vertical  asymmetric  inlets  in  order  to  induce  entrainment  flow  and  generate  a 
swirling  motion.  The  model  predictions  gave  quantitatively  similar  results  to  that 
observed  in  experiments;  however,  the  simulations  did  not  include  a combustion  model. 

Among  the  studies  directed  specifically  toward  fire  whirls,  there  also  exists  research 
in  the  areas  of  swirling  [9,  10]  and  whirling  [11]  non-premixed  flames.  According  to 
Gabler  et  al.  [11],  whirling  flames  differ  from  swirling  flames  in  that  swirl  combines 
both  axial  and  tangential  motions  of  a flow.  Whirling  flames  are  described  as  a rotation 
produced  purely  by  tangential  motion.  The  studies  of  swirling  flames  include  details 
on  the  effects  of  varying  swirl  and  jet  momentum.  The  fire  whirl  experiment  [1]  was 
extended  by  Chigier  et  al.  [9]  whereby  the  liquid  pool  was  replaced  with  a turbulent 
jet  diffusion  flame.  The  experiments  yielded  increased  flame  lengths  of  the  fuel  jet  with 
either  increasing  rotation,  increasing  flow  rates,  or  both.  Tangirala  and  Driscoll  [10] 
proceeded  with  studies  of  a fuel  jet  surrounded  by  coaxial  swirling  air.  It  was  reported 
that  increasing  swirl  decreased  the  jet-like  flame  to  a short  intense  flame,  with  improved 
stability,  rapid  mixing,  and  increased  temperatures.  Increasing  fuel  jet  momentum 
produced  a more  jet-like  flame.  At  this  point,  it  is  important  to  recognize  that  jet 
flames  [9,  10,  11]  behave  differently  depending  on  the  application  of  a rotating  flowfield. 
It  would  seem  that  flame  lengths  tend  to  decrease  when  rotation  is  imparted  to  the  jet, 
but  increase  when  rotation  is  imparted  to  the  environment. 

To  date,  fire  whirls  have  not  been  adequately  studied  either  theoretically  or 
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experimentally,  so  the  impact  of  swirl  on  the  fluid  dynamics  and  combustion  process 
is  not  fully  understood.  Reproducing  the  characteristics  of  a fire  whirl  for  any  study, 
whether  it  be  experimental  or  numerical,  is  an  issue  that  is  addressed  in  this  paper.  The 
swirling  motion  of  a fire  has  been  initiated  within  a rotating  cylindrical  screen  [1]  and 
vented  walls  [7,  8].  However,  if  one  is  to  understand  the  relationship  between  buoyancy 
and  angular  momentum,  it  is  necessary  to  identify  the  magnitude  of  each  attribute. 
We  therefore  approach  the  problem  of  simulating  a fire  whirl  by  imposing  circulation 
into  the  buoyant  plume  scenario.  In  this  manner  we  are  able  to  vary  parameters  which 
provides  a controlled  analysis  of  the  fundamental  dynamics  of  swirling  fire  plumes. 

The  phenomenon  of  a swirling  fire  can  be  described  physically  with  critical 
parameters  that  control  the  dynamic  mechanisms.  To  begin,  the  overall  heat  release 
rate  determines  the  strength  of  the  fire.  In  the  absence  of  any  swirl,  the  fire  contains 
a natural  length  scale  that  is  a characteristic  flame  height  based  on  the  overall  heat 
release  rate  and  ambient  air  properties.  The  length  scale  together  with  the  gravitational 
acceleration  are  used  to  define  scales  for  velocity  and  time.  The  entire  problem  can  be 
made  dimensionless  such  that  the  isolated  plume  depends  only  on  a Reynolds  number 
characterized  by  an  effective  kinematic  viscosity.  The  modifications  introduced  by 
imposing  circulation  can  now  be  related  to  the  defined  characteristic  scales.  It  is  the  ratio 
of  the  imposed  (swirl)  circulation  to  the  buoyancy  induced  circulation  that  characterizes 
the  relative  importance  of  these  mechanisms. 

The  experiments  described  above  indicate  three  possible  kinds  of  scenarios.  First, 
the  swirl  induced  circulation  is  held  fixed,  and  the  fire  strength  is  allowed  to  change 
[1].  Second,  the  fire  strength  is  fixed  and  the  geometry  determines  the  strength  of  the 
induced  swirl  [7].  Finally,  it  is  possible  to  control  both  the  swirl  induced  circulation  and 
the  fire  strength  using  a rotating  screen  to  generate  the  swirl  and  a burner  (as  opposed 
to  a pool  fire)  to  generate  buoyancy  [9].  This  latter  scenario  can  be  considered  as  the 
“baseline”  fire  whirl,  whose  dynamics  is  the  subject  of  the  present  paper.  Figure  1 is  a 
schematic  of  the  configuration.  It  is  hoped  that  these  constraints  will  be  systematically 
relaxed  in  future  papers  using  the  methods  described  below. 

The  numerical  investigations  of  swirling  fire  plumes  are  used  to  analyze  how  swirl 
alters  the  plume  dynamics  and  combustion.  A methodology  has  been  developed  to  model 
fires  with  the  use  of  efficient  flow  solving  techniques  and  high  spatial  resolution  [12,  13]. 
An  approximate  form  of  the  Navier-Stokes  equations  is  solved  to  calculate  the  mixing 
and  transport  of  combustion  products.  Large  scale  eddies  are  directly  simulated  and 
sub-grid  scale  motion  is  represented  with  a Smagorinsky  model.  The  fire  is  prescribed 
in  a manner  consistent  with  a mixture-fraction  based  approach  to  combustion  [14,  15]. 
The  models  and  computational  methodology  have  reproduced  mean  temperature  and 
buoyant  velocity  correlations  for  large  fire  plumes  in  the  absence  of  circulation  [16,  17]. 

The  paper  is  laid  out  as  follows.  Section  2 begins  with  a description  of  the  problem, 
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Figure  1.  Schematic  of  a fire  plume  in  ambient  conditions  under  an  imposed 
circulation  To  (e.g.,  a rotating  cylindrical  screen).  The  fire  bed  is  of  heat  release 
rate  q with  height  zq  and  flame  height  D* . 

and  current  large  eddy  simulation  (LES)  methodology  and  combustion  model.  Section  3 
is  a presentation  and  discussion  of  the  LES  results  for  a fire  plume  in  ambient  conditions 
with  and  without  swirl.  The  discussion  begins  with  a comparison  of  non-swirling  fires 
and  well-established  plume  correlations.  Both  non-swirling  and  swirling  plumes  are 
then  contrasted  to  highlight  basic  physical  differences.  An  important  contribution  of 
this  paper  is  the  presentation  of  temperature  fields,  swirl  velocity  characteristics,  and 
heat  release  rate  trends  with  increasing  circulation.  The  summary  will  address  the 
advantages  and  limitations  of  the  simulations  to  establish  goals  of  future  work. 

2.  Methodology 

2.1.  Hydrodynamic  model 

The  fire  plume  is  a three-dimensional,  transient,  buoyant  flow  that  can  be  modeled  by 
the  motion  of  a thermally  expandable  ideal  gas  [12].  The  Navier-Stokes  equations  are 
solved  for  such  a fluid  driven  by  a prescribed  heat  source,  where 


dp 

- + V.pu  = ° 

(1) 

^ + V ■ pYiU  = V ■ pDiVYi  + W 
dt 

(2) 

(du  \ 

pi  — + u-Vul  + Vp  - pg  = V ■ t 

(3) 
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~ + V . phu  _ ^ = g'"  + V • (fcVT  + £ AAVKi  + ?«)  (4) 

p0(t)  = pTH^Yt/Mi)  (5) 

The  fluid  variables  are  density  p,  velocity  u,  mass  fraction  of  zth  species  M,  mass 
diffusivity  of  2th  species  Dr , production  rate  of  2th  species  ITj,  pressure  p,  gravity  g, 
viscous  stress  tensor  r,  enthalpy  h.,  volumetric  heat  release  rate  qm,  thermal  conductivity 
k,  temperature  T,  radiant  energy  flux  qR , universal  gas  constant  7Z , and  molecular  weight 
of  2th  species  M*. 

The  pressure  is  decomposed  into  three  components,  a background  (average)  pressure 
Po(t),  a hydrostatic  contribution  — p0gz , and  a perturbation  to  the  hydrostatic  p(x,t). 
The  subscript  0 is  a reference  value  of  density,  and  z is  the  vertical  coordinate.  Equations 
(4)  and  (5)  are  modified  using  the  background  pressure  po  which  depends  on  time  t 
only.  In  this  manner,  high-frequency  acoustic  oscillations  are  eliminated  while  large 
temperature  and  density  variations  typically  found  in  fires  are  retained.  The  resulting 
equations  are  referred  to  as  weakly  compressible  and  are  valid  for  low  Mach  number 
flows  [12]. 

A further  assumption  is  that  the  constant  pressure  specific  heat  of  the  2th  species 
cPii  is  independent  of  temperature.  Therefore,  the  enthalpy  is  written  as 

h = '£ihiYi  = TY,Cp,iYi 

Combining  the  conservation  equations  for  mass  (1),  energy  (4)  and  state  (5)  results  in 
an  expression  for  the  divergence  of  the  flow  V • u 

Po(t)V  ■ u + [q"'  + V • [kVT  + Y,ph,D,VY,  + <?*)]  (6) 

The  term  V-u  is  an  important  quantity  that  is  later  exploited  in  the  numerical  solution. 
Equation  (6)  is  integrated  over  the  entire  domain  for  pressure  p0  and  is  used  as  a 
consistency  condition. 

The  momentum  equation  (3)  is  simplified  by  subtracting  off  the  hydrostatic  pressure 
gradient,  and  then  dividing  by  the  density  to  obtain 

du  1„,  ,9  1 1 ,,  , _ . 

"XT  + o Vlul  - U X + — Vp  =-((p-  p0)g  + V • t) 
at  2 po  p 

where  u>  is  the  vorticity.  The  head  H is  introduced,  where 

V%  = t:V|u|2  + — Vp 

2 Po 

The  final  form  of  the  momentum  equation  is 

O ^ 1 

"UXW  + vn  = - ((p  - Po)g  + v • r) 

ot  p 


(7) 
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The  quantity  1 / p0  implies  that  the  vorticity  generation  induced  by  buoyancy  is  much 
more  important  than  that  due  to  baroclinic  effects,  i.e.  the  non-alignment  of  the  pressure 
and  density  gradients. 

The  equations  presented  thus  far  can  be  solved  directly  to  obtain  a solution  to  the 
problem  of  a buoyant  fire  plume.  The  difficulty  arises  in  that  the  length  and  time  scales 
associated  with  the  fluid  dynamics  and  combustion  vary  over  orders  of  magnitude.  In 
order  to  calculate  the  dynamics  of  the  problem,  models  for  r and  q"'  are  necessary. 
The  goal  is  to  capture  mixing  on  scales  where  the  eddy  viscosity  is  effective  which 
generally  requires  a very  fine  grid  to  resolve  the  scales.  An  alternative  approach  is  to 
use  an  elaborate  sub- grid  scale  description.  One  representation  of  the  dynamic  viscosity 
is  based  on  the  analysis  of  Smagorinsky  [18]  where  the  sub-grid  scale  Reynolds  stress 
tensor  is  given  by 


Aj  P 


+ 


9Uj 

dxi 


P — P {Cs  A)2 15| 


where  Cs  is  an  empirical  constant,  A is  a length  on  the  order  of  the  grid  cell  size,  and 


n (dw' 
+ 2{dz 


+ 


du  dv 
dy  + dx 


( du  dw  \ 

+ + dx  J 


( dv  dw' 

1 + 


There  have  been  numerous  refinements  of  the  original  Smagorinsky  model  but  it  is 
difficult  to  assess  the  improvements  offered  by  these  newer  schemes.  There  are  two 
reasons  for  this.  First,  the  structure  of  the  fire  plume  is  dominated  by  the  large  scale 
resolvable  eddies  so  that  even  a constant  eddy  viscosity  gives  results  almost  identical 
with  those  obtained  with  the  Smagorinsky  scheme  [19].  Second,  the  lack  of  precision 
in  most  large  scale  fire  data  makes  it  difficult  to  sort  out  the  subtleties  associated  with 
these  models.  For  the  time  being,  the  Smagorinsky  model  with  Cs  — 0.14  produces 
satisfactory  results  for  most  large  scale  applications  where  boundary  layers  are  not 
important. 

The  pressure  perturbation  is  solved  by  taking  the  divergence  of  the  momentum 
equation  which  results  in  an  elliptic  partial  differential  equation  for  H 

9 <9(V  • u) 

V2H  ^ - V • F 

dt 

where  the  convective  and  diffusive  terms  have  been  incorporated  in  the  term  F. 


2.2.  Combustion  model 

For  the  large  eddy  simulations,  a sub-grid  thermal  element  model  (TEM)  is  formulated 
to  represent  the  fire.  A large  number  of  Lagrangian  particles  are  introduced  into  the 
plume,  releasing  heat  as  they  are  convected  by  the  thermally  induced  motion  [13,  14]. 
The  combustion  and  hydrodynamics  are  coupled  since  the  fluid  motion  determines  where 
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the  heat  is  released,  while  the  heat  release  determines  the  motion.  The  concept  is 
based  on  mixture-fraction  theory  for  the  transport  of  a conserved  scalar  that  describes 
a reacting  species.  The  TEM  requires  that  the  burning  rate  be  prescribed  as  an  input 
into  the  calculation.  It  is  intended  for  applications  where  small  scale  mixing  and  diffusive 
processes  that  control  non-premixed  combustion  cannot  be  resolved  on  the  grid  used  to 
perform  the  simulation. 

The  overall  heat  release  rate  q from  the  fire  is  discretized  as  thermal  elements  that 
represent  pyrolyzed  fuel.  At  a specified  surface,  such  as  the  fuel  bed,  thermal  elements 
are  ejected  at  a rate  of  h"  particles  per  unit  time  per  unit  area  with  a small  normal 
velocity  into  the  flow  domain.  The  heat  release  rate  of  a single  thermal  element  is  given 
by 

q"  1 

Qp,j  — - n , \t  t0  < tb) 

Tl  Zb 

where  q"  is  the  heat  release  rate  per  unit  area  of  the  fuel  bed  and  tb  is  the  burn-out  time 
of  the  thermal  element.  The  burn-out  time  is  obtained  from  the  plume  correlations  of 
Baum  and  McCaffrey  [17].  It  is  assumed  that  the  thermal  element  burns  out  somewhere 
in  the  intermittent  region  of  the  plume,  1.32  D*  < z < 3.30  .D*,  where  z is  the  height 
above  the  fire  bed,  D*  = (q/ {poCpT0y/g))2t5  is  the  characteristic  diameter  of  the  fire,  and 
q is  the  total  heat  release  rate  of  the  fire.  The  burn-out  time  falls  somewhere  between 
1.05 yj D* / g < tb  < 1.8 6^D* /g  and  is  usually  a few  tenths  of  a second.  The  heat  release 
term  in  (4)  is  the  summation  of  convective  heat  release  rates  of  the  individual  thermal 
elements  in  a grid  cell  of  volume  5xSySz. 


Y Qp,j 

SxSySz 


For  the  inclusion  of  oxygen  transport  in  the  calculation,  the  burn-out  time  of  any 
thermal  element  will  vary  based  on  the  concentration  of  oxygen  in  the  surrounding  gas. 
Oxygen  is  consumed  in  any  given  control  volume  based  on  the  amount  of  heat  generated 
in  the  control  volume.  The  source  term  in  the  oxygen  transport  equation  (2)  becomes 


W02  = 


AH, 


o2 


where  A Hq2  is  the  amount  of  heat  liberated  per  unit  mass  of  oxygen  consumed  (usually 
about  13,100  kJ/kg  O2).  When  the  oxygen  mass  fraction  Yq2  falls  to  a certain  prescribed 
lower  limit,  combustion  is  assumed  to  stop,  and  the  unburned  fuel  associated  with  the 
thermal  elements  remains  unburned  until  more  oxygen  is  available. 


2.3.  Numerical  formulation 

The  equations  numerically  solved  are  the  conservation  of  mass  (1),  species  (2), 
momentum  (3),  and  the  Poisson  equation  for  the  total  pressure  which  is  formulated  from 
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the  energy  equation  (4)  and  state  equation  (5).  The  combustion  is  simulated  via  the 
thermal  element  model  as  described  in  section  2.2.  The  spatial  derivatives  are  discretized 
with  second-order  central  differencing,  and  a second-order  explicit  Runge-Kutta  scheme 
is  used  to  advance  the  velocity  and  temperature  fields.  The  linear  algebraic  system 
arising  from  the  discretization  of  the  Poisson  equation  (8)  has  constant  coefficients 
and  can  be  solved  to  machine  accuracy  by  direct  (non-iterative)  methods  that  utilize 
fast  Fourier  transforms  and  block  tri-diagonal  solvers.  The  grid  is  rectangular  with 
uniform  spacing  in  the  horizontal  directions  and  stretching  in  the  vertical  direction. 
The  discretization  for  a computational  cell  is  based  on  staggered-gridding  techniques. 

No-flux  boundary  conditions  are  specified  by  asserting  that 


&H 

dn 


= -Fn 


at  solid  walls,  where  Fn  is  the  normal  component  of  F at  the  wall.  Therefore,  the  normal 
component  of  velocity  at  the  wall  does  not  change  with  time,  and  indeed  remains  zero 
assuming  the  flow  velocity  is  initially  zero.  At  open  external  boundaries  it  is  assumed 
that  the  perturbation  pressure  is  zero.  Thermal  elements  are  ejected  from  the  burner 
surface  and  burn  according  to  the  heat  release  rate  per  unit  area  q"  which  is  specified 
as  an  input.  An  adiabatic  boundary  condition  is  used  so  that  there  is  no  temperature 
gradient  normal  to  the  burner  surface.  Refer  to  McGrattan  et  al.  [13]  for  further  details. 


3.  Large  eddy  simulations 

The  simulations  of  a whirling  fire  are  based  on  the  schematic  shown  in  figure  1.  The 
domain  is  a rectangular  region  with  a square  base  0.50  m per  side  and  a 1.0  m height. 
A circular  burner  is  centered  at  the  base  of  the  domain  with  diameter  d = 0.0863  m and 
height  z0  = 0.05  m.  The  origin  of  the  Cartesian  coordinate  system  is  centered  at  the 
base.  The  burner  specifications  represent  a fuel  source  of  acetone  yielding  a 14.4  kW 
fire.  The  oxygen  transport  model  was  used  with  a prescribed  mass  fraction  burnout  of 
0.15. 

Analogous  to  the  experiments  of  Emmons  and  Ying  [1],  swirl  can  be  induced  by 
simulating  an  imposed  circulation  r0  on  the  fire  plume  at  a prescribed  radial  distance  r0 
from  the  nominal  centerline.  Three  dimensionless  parameters  govern  the  flow  generated 
by  this  model: 

(i)  the  heat-release  rate  q made  dimensionless  using  thermal  properties  of  the  ideal  gas 

(ii)  the  circulation  made  dimensionless  using  a characteristic  flame  diameter  D*  [17] 
and  buoyant  velocity  \/gD*,  where  Q/a  = ( r0/27rrl)/a 

(iii)  a Reynolds  number  based  on  the  buoyant  velocity  and  flame  diameter. 
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The  parameter  a is  the  reciprocal  of  a time  scale  which  characterizes  the  flame  strain 
due  to  buoyancy  a = yj g/ D* . Therefore,  the  parameter  Q/a  provides  a measure  of 
the  effects  of  swirl-induced  motion  to  buoyancy-driven  motion.  The  Reynolds  number 
is  defined  using  characteristic  parameters  most  relevant  to  the  described  problem. 
However,  its  definition  should  not  be  confused  with  the  dynamic  Reynolds  number  of  the 
large  scale  simulations  (refer  to  section  2.1).  Additional  dimensionless  parameters  are 
the  ratio  of  length  scales  r0/D*  and  d/D *;  note  that  for  these  simulations  the  parameters 
are  fixed  and  of  order  unity. 

Results  are  presented  for  the  large  eddy  simulations  using  a rectangular  grid. 
Numerical  tests  were  performed  to  determine  the  grid  size  necessary  to  resolve  the 
flowfield  with  the  best  accuracy.  Five  mesh  sizes  ranging  from  323  cells  to  1283  cells 
were  used  in  the  analysis.  Based  on  Richardson’s  extrapolation , simulations  calculated 
on  a 723  grid  (consisting  of  373,248  cells)  yields  a relative  error  of  about  4%.  An  average 
simulation  takes  approximately  18  ps  per  time  step  per  cell  on  an  SGI  Octane/SI  R10000 
195MHz  processor, J and  requires  about  12000  time  steps  for  20  s of  simulation  time. 

The  LES  model  was  found  to  reproduce  mean  temperature  and  buoyant  velocity 
correlations  for  large  fire  plumes  in  the  absence  of  circulation  over  most  of  the  simulated 
domain.  Figures  2(a-b)  compare  the  simulations  with  the  McCaffrey  correlations  for 
the  time-averaged  centerline  temperature  and  vertical  velocity,  respectively.  Baum  and 
McCaffrey  [17]  define  three  major  zones  for  a typical  fire  plume.  The  zones  are  referred 
to  as  1.)  the  continuous  (visible)  flame  for  0 < (z  — z0)/D*  < 1.32,  2.)  the  intermittent 
region  for  1.32  < (z  — z0)/D*  < 3.30,  and  3.)  the  plume  region  for  3.30  < (z  — z0)/D*. 
Time-averaged  radial  temperature  and  velocity  profiles  have  been  shown  previously 
by  Baum,  et  al  [19].  The  good  correspondence  between  the  correlation  data  and  the 
simulations  is  very  encouraging,  and  provides  confidence  for  the  LES  formulation. 


Figure  2.  Comparison  of  the  McCaffrey  correlations  [17]  with  LES  for  the  time- 
averaged  centerline  a)  temperature,  and  b)  vertical  velocity. 
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However,  it  is  worthwhile  to  note  that  near  the  base  of  the  fire  the  correspondence 
between  the  temperature  correlations  and  LES  deviates.  The  discrepancy  is  not  a failure 
of  the  correlations  but  is  a limitation  of  the  current  model.  Neither  the  correlations  nor 
the  numerical  model  were  constructed  to  apply  near  the  base  of  a fire,  but  rather  serve  as 
a global  description  of  flow  dynamics.  There  is  presently  insufficient  concensus  between 
experiments  to  properly  guide  the  formulation  of  a model  to  represent  the  entrainment 
region  close  to  the  fuel  source  [20].  The  reason  for  this  is  two-fold:  1.)  the  region 
near  the  base  of  the  fire  is  sensitive  to  the  fuel  geometry  and  can  drastically  change 
depending  on  the  boundary  conditions,  experimental  apparatus,  etc.,  and  2.)  there  is 
great  uncertainty  in  the  measurements  of  entrainment  rates  near  the  fuel  source  due  to 
localized  effects. 

The  fire  plume  was  simulated  with  an  imposed  circulation  at  a distance  r0  = 0.2  m 
measured  from  the  nominal  plume  center  and  along  the  entire  height  of  the  domain.  A 
comparison  of  the  instantaneous  temperature  fields  for  a non-swirling  (fi/c*  = 0.0)  and 
swirling  fire  {Q./a  — 0.536)  is  shown  in  figures  3(a-b),  respectively.  The  non-swirling 
flame  in  figure  3(a)  shows  the  onset  of  characteristic  puffing  associated  with  a typical 
fire  [19,  21]  as  two  hot  regions  near  z/D*  — 0.5  and  1.0.  However,  the  swirling  fire 
(figure  3(b))  tends  to  have  a continuous  flame  up  to  z/D*  — 2.3.  Another  difference  is 
that  the  temperature  contours  shown  in  figure  3(b)  define  the  continuous  flame  as  being 
more  tapered  and  constricted  near  the  centerline  than  the  non-swirling  flame.  Contours 
of  heat  release  rate  per  unit  volume  q'"  are  shown  in  figure  4(a-b)  for  Q/a  — 0.0  and 
0.536.  Note  that  the  swirling  fire  has  more  regions  of  high  heat  release  rate  than  the 
non-swirling  fire. 

A three-dimensional  representation  of  the  two  cases  is  shown  in  figures  5(a-b).  The 
figures  are  isotherms  for  three  temperature  levels  corresponding  to  the  different  plume 
regimes  defined  previously.  Sample  Lagrangian  trajectories  are  shown  as  ribbons  which 
follow  the  path  of  a streamline.  The  trajectory  of  figure  5(a)  shows  how  a particle 
entrained  in  the  non-swirling  fire  is  drawn  towards  the  center  of  the  plume  and  then 
merely  travels  straight-up.  A particle  entrained  in  a swirling  fire  is  also  drawn  towards 
the  center  of  the  plume  but  wraps  around  the  plume  a few  times  near  the  base  of  the 
fire  before  accelerating  vertically.  The  general  fluid  motion  of  a swirling  fire  is  unstable, 
as  the  fire  typically  oscillates  or  precesses  within  the  boundaries  of  the  fire  bed.  Bearing 
this  in  mind,  time-averaged  data  is  now  presented. 

The  time-averaged  temperature  fields  shown  in  figures  6(a-b)  are  not  as  elucidating 
as  the  instantaneous  fire  images  (3(a-b)).  The  temperature  contours  of  the  non-swirling 
fire  (3(a))  are  more  symmetric  than  the  swirling  fire  (3(b)).  It  appears  that  time- 
averaging smears  out  the  precessing  motion  of  the  swirling  flows.  However,  to  draw 
comparisons  with  experimental  work,  the  remainder  of  this  paper  will  consider  results 
for  the  mean  flowfields. 
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Figure  3.  Instantaneous  temperature  field  at  the  plume  centerplane  for  a)  non-swirling 
fire  f }/a  = 0.0,  and  b)  whirling  fire  Cl/a  = 0.536. 
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Figure  4.  Instantaneous  heat  release  rate  per  unit  volume  at  the  plume  centerplane 


for  a)  non-swirling  fire  f l/a  = 0.0,  and  b)  whirling  fire  Cl /a  = 0.536. 


Figure  5.  Simulations  of  a plume  for  a)  non-swirling  fire  Cl/a  = 0.0,  b)  whirling  fire 
Cl/a  = 0.536.  Isotherms  are  shown  for  three  major  zones:  the  continuous  flame  (dark 
gray),  the  intermittent  region  (medium  gray),  and  the  plume  region  (light  gray).  The 
black  ribbons  are  sample  Lagrangian  trajectories  of  the  flowfields. 


Figure  6.  Time-averaged  temperature  fields  at  the  plume  centerplane  for  a)  non- 
swirling fire  Cl /a  = 0.0,  and  b)  whirling  fire  Cl /a  = 0.536. 
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Figure  7.  Time-averaged  temperature  fields  at  the  plume  centerplane  for  varying 
circulation  f 1/a. 


Figure  7 shows  six  cases  of  the  time-averaged  temperature  fields  for  values  of 
circulation  from  Q/a  = 0.0  to  Q/a  ~ 0.8.  At  first,  increasing  the  circulation  results 
in  a larger  region  of  high  temperatures  near  the  fire  bed.  The  simulations  then  predict 
that  further  increasing  the  circulation  above  Q/a  > .268  reduces  the  region  of  high 
temperatures.  However,  the  isotherms  near  the  central  region  of  the  plume  show  that 
the  flame  has  narrowed  and  stretched  vertically  for  Q/a  > 0.536.  The  flame-lengthening 
trends  discussed  by  Chigier  et  al.  [9}  are  qualitatively  similar  to  that  reported  here. 
Furthermore,  there  is  very  good  agreement  when  cases  Q/a  = 0.0  and  0.804  (figure  7) 
are  compared  with  the  temperature  distributions  of  the  propane  diffusion  flame  [9],  for 
both  flame  shape  and  temperature  ranges. 

Another  way  to  interpret  the  varying  effects  of  circulation  is  shown  as  contours  of 
time-averaged  heat  release  rate  per  unit  volume  in  figure  8.  Increasing  the  circulation 
reduces  the  overall  volume  of  the  fire  plume  and  increases  the  region  of  high  heat  release 
rates  within  the  central  region  of  the  flame.  The  heat  release  rate  contours  near  the 
base  for  Q/a  > 0.50  have  two  peaks  and  a small  dip  near  the  centerline,  which  seems  to 
be  a consequence  of  time-averaging  an  oscillating  flow.  The  experimental  observations 
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Figure  8.  Time-averaged  heat  release  rate  per  unit  volume  at  the  plume  centerplane 
for  varying  circulation  Q/a. 


of  Emmons  and  Ying  [1]  indicate  that  the  heat  release  rate  is  significantly  increased  by 
the  swirling  motion.  Although  the  numerical  methodology  described  in  section  2 does 
not  model  varying  heat  release  rate,  it  is  very  encouraging  that  the  simulated  results 
show  larger  regions  of  higher  heat  release  rates. 

Time-averaged  temperature  profiles  for  the  six  cases  of  circulation  are  shown  in 
figure  9 at  (z  — z0)/D*  — 1.5.  Overall,  the  temperature  profiles  have  nearly  the  same 
distribution  along  x at  the  centerplane.  A closer  examination  of  the  temperatures  near 
the  centerline  {x  = 0.0)  reveals  an  interesting  pattern.  Initially,  the  mean  temperature  at 
the  centerline  reduces  with  increasing  circulation  up  to  Q/a  = 0.536.  Further  increasing 
circulation  then  results  in  an  increase  of  the  mean  temperature  at  the  centerline.  The 
decreasing/increasing  temperature  trends  are  also  reflected  in  figure  7 of  time-averaged 
isotherms. 

At  this  point  we  explore  the  possibility  of  defining  a characteristic  flame  height 
for  the  six  cases  presented.  According  to  analyses  of  Orloff  and  deRis  [22],  the  size 
and  shape  of  pool  fires  can  be  constructed  through  Froude  modeling.  Their  model  was 
substantiated  with  experimental  data  from  pool  fires  of  various  fuels.  The  luminous 
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Figure  9.  Mean  temperature  profiles  at  (z  — zq)/ D*  = 1.5  for  varying  circulation. 
Figure  10.  Mean  centerline  flame  height  (based  on  925° C)  versus  circulation. 
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flame  volume  was  defined  by  an  effective  flame  temperature  of  1200  K (independent  of 
fuel)  for  which  the  Froude  modeling  suggested  an  average  volumetric  heat  release  rate 
of  1200  kW/m3.  Returning  to  figure  7,  the  flame  height  can  be  visualized  as  the  contour 
marking  925°C.  Likewise,  figure  8 deliniates  the  flame  shape  between  contours  1000-1500 
kW/m3.  A characteristic  mean  flame  height  at  the  centerline  based  on  925°C  is  shown 
in  figure  10  to  re-emphasize  the  decreasing/increasing  flame  height  trends  previously 
discussed. 

The  numerical  simulations  can  be  compared  to  the  experimental  findings  of  Emmons 
and  Ying  [1].  Their  study  of  a fire  whirl  included  temperature  profiles  for  different  values 
of  circulation.  The  experimental  trends  indicate  that  with  increasing  circulation,  the 
maximum  mean  temperature  peaked  off-center.  Figures  ll(a-b)  are  mean  temperature 
profiles  for  Q/a  = 0.0  and  0.536.  The  maximum  mean  temperature  for  the  non-swirling 
fire  always  occurs  at  the  plume  nominal  centerline  and  decreases  with  height.  In  the  case 
of  the  swirling  fire,  the  maximum  temperatures  also  occurred  at  the  nominal  centerline. 
Further  data  analyses  indicated  that  depending  on  the  number  of  realizations  captured 
in  one  complete  revolution  of  the  fire,  the  time-averaged  data  may  yield  asymmetric 
features  of  the  flow.  It  is  not  clear  if  the  off-centerline  peak  temperatures  as  reported 
in  [1]  are  an  artifact  of  the  oscillating  motion  of  the  fire  or  due  to  the  location  of  the 
flame  sheet. 

The  swirl  velocity  profiles  provide  further  information  regarding  the  fluid  motion. 
The  time- averaged  profiles  are  shown  in  figure  12  for  varying  circulation  at  position 
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(. z—zq)/D * = 1.5.  In  the  absence  of  circulation,  the  mean  swirl  velocity  is  approximately 
zero  with  some  fluctuation  near  the  plume  centerline  and  closest  to  the  fire  bed.  Once 
circulation  is  introduced,  the  magnitude  of  the  swirl  velocity  increases  with  increasing 
circulation.  Peak  values  of  swirl  velocity  occur  at  x/D*  ~ 0.4  which  corresponds  to 
the  edge  of  the  burner  and  may  be  indicative  of  high  heat  release  rates  near  the  edge 
of  the  flame.  Beyond  x/D*  = 1.1,  the  position  at  which  circulation  is  imposed,  the 
swirl  velocity  decays  to  zero.  Figure  13  is  a plot  of  the  mean  swirl  velocity  versus 
x/D*  at  various  heights  for  Q/a  = 0.536.  The  curves  indicate  that  the  peak  swirl 
velocity  decreases  with  increasing  height  z.  As  the  plume  spreads  radially,  the  swirl 
velocities  tend  to  mix  outward  and  reduce  the  swirl  velocities,  suggesting  that  buoyancy 
increasingly  dominates  the  dynamics  farther  up  the  plume. 

4.  Concluding  remarks 

Numerical  simulations  of  fire  plumes  were  conducted  in  order  to  ascertain  the  dynamic 
attributes  of  swirl  and  buoyancy  in  whirling  fires.  The  LES  methodology  for  simulating 
fires  is  based  on  direct  solution  of  the  large  scales  and  sub-grid  scale  approximations 
for  the  dynamic  viscosity  and  heat  release.  The  simulations  for  non-swirling  fires  are  in 
very  good  agreement  with  mean  temperature  and  buoyant  velocity  correlations  for  large 
fires.  The  approach  taken  for  swirling  fires  was  to  impose  circulation  while  maintaining 
a constant  heat  release  rate. 

It  was  shown  that  increasing  circulation  has  a strong  influence  on  the  shape  of  the 
fire  plume.  The  whirling  fires  increased  in  length  and  constricted  radially,  consistent 
with  previously  published  literature.  The  temperature  trends  along  the  centerline 
decreased  then  increased  in  values  for  increasing  circulation.  Despite  the  fixed  heat 
release,  increasing  circulation  increased  the  region  of  maximum  heat  release  rates  while 
reducing  the  overall  volume  of  the  fire.  The  swirl  velocities  also  increased  with  increasing 
swirl.  In  addition,  it  was  determined  that  farther  up  the  plume  buoyancy  effects  are 
more  dominant  than  the  swirl  velocities. 

The  present  study  carefully  explored  the  effects  of  circulation  on  a fire  plume.  The 
scenario  of  a fixed  fire  strength  and  variable  circulation  was  the  first  numerical  study 
to  systematically  analyze  the  relationship  between  swirl  and  buoyancy.  A set  of  well- 
defined  experiments  to  match  the  current  numerical  approach  will  be  beneficial  for 
a more  quantitative  comparison  and  validation.  The  next  step  will  be  to  relax  the 
condition  of  a fixed  heat  release  to  one  that  varies  according  to  the  radiation  heat 
feedback  to  the  surface  of  the  fire.  This,  in  itself,  changes  the  current  problem  of  an 
imposed  fuel  flow  to  a pool  fire. 

t The  computer  is  identified  to  specify  the  computational  timings  and  is  not  intended  to  imply 
recommendation  nor  endorsement  by  NIST. 
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(b) 

Figure  11.  Time-averaged  temperature  profiles  for  a)  non-swirling  fire  Q/a  — 0.0 
and  b)  whirling  fire  Cl/a  = 0.536,  for  varying  height  (z  — zq)/D*. 


Figure  12 


Figure  13 

Figure  12.  Mean  swirl  velocity  profiles  at  (z  - zq)/D*  — 1.5  for  varying  circulation. 
Figure  13.  Mean  swirl  velocity  profiles  of  whirling  fire  f 1/a  = 0.536  for  varying  height. 
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